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1. Preamble 

These lecture notes present the computation of the full system of Boltzmann 
equations describing the evolution of the photon, baryon and cold dark matter flu- 
ids up to second order in perturbation theory, as recently studied in Refs. [1,2]. 
These equations allow to follow the time evolution of the cosmic microwave 
background anisotropics at all angular scales from the early epoch, when the 
cosmological perturbations were generated, to the present, through the recombi- 
nation era. The inclusion of second-order contributions is mandatory when one 
is interested in studying possible deviations from Gaussianity of cosmological 
perturbations, either of primordial (e.g. inflationary) origin or due to their sub- 
sequent evolution. Most of the emphasis in these lectures notes will be given 
to the derivation of the relevant equations for the study of cosmic microwave 
background anisotropies and to their analytical solutions. 



2. Introduction 

Cosmic Microwave Background (CMB) anisotropies play a special role in cos- 
mology, as they allow an accurate determination of cosmological parameters and 
may provide a unique probe of the physics of the early universe and in particular 
of the processes that gave origin to the primordial perturbations. 

Cosmological inflation [3] is nowadays considered the dominant paradigm for 
the generation of the initial seeds for structure formation. In the inflationary pic- 
ture, the primordial cosmological perturbations are created from quantum fluctu- 
ations "redshifted" out of the horizon during an early period of accelerated expan- 
sion of the universe, where they remain "frozen". They are observable through 
CMB temperature anisotropies (and polarization) and the large-scale clustering 
properties of the matter distribution in the Universe. 

This picture has recently received further spectacular confirmations from the 
results of the Wilkinson Microwave Anisotropy Probe (WMAP) three year set 
of data [4]. Since the observed cosmological perturbations are of the order of 
10 -5 , one might think that first-order perturbation theory will be adequate for all 
comparisons with observations. This might not be the case, though. Present [4] 
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and future experiments [5] may be sensitive to the non-linearities of the cosmo- 
logical perturbations at the level of second- or higher-order perturbation theory. 
The detection of these non-linearities through the non-Gaussianity (NG) in the 
CMB [6] has become one of the primary experimental targets. 

There is one fundamental reason why a positive detection of NG is so relevant: 
it might help in discriminating among the various mechanisms for the generation 
of the cosmological perturbations. Indeed, various models of inflation, firmly 
rooted in modern particle physics theory, predict a significant amount of primor- 
dial NG generated either during or immediately after inflation when the comov- 
ing curvature perturbation becomes constant on super-horizon scales [6]. While 
single-field [7] and two(multi)-field [8] models of inflation predict a tiny level of 
NG, "curvaton"-type models, in which a significant contribution to the curvature 
perturbation is generated after the end of slow-roll inflation by the perturbation 
in a field which has a negligible effect on inflation, may predict a high level of 
NG [9]. Alternatives to the curvaton model are models where a curvature per- 
turbation mode is generated by an inhomogeneity in the decay rate [10, 1 1], the 
mass [12] or the interaction rate [13] of the particles responsible for the reheating 
after inflation. Other opportunities for generating the curvature perturbations oc- 
cur at the end of inflation [14], during preheating [15], and at a phase-transition 
producing cosmic strings [16]. 

Statistics like the bispectrum and the trispectrum of the CMB can then be 
used to assess the level of primordial NG on various cosmological scales and to 
discriminate it from the one induced by secondary anisotropics and systematic 
effects [6, 17-19]. A positive detection of a primordial NG in the CMB at some 
level might therefore confirm and/or rule out a whole class of mechanisms by 
which the cosmological perturbations have been generated. 

Despite the importance of NG in CMB anisotropics, little effort has ben made 
so far to provide accurate theoretical predictions of it. On the contrary, the vast 
majority of the literature has been devoted to the computation of the bispectrum 
of either the comovig curvature perturbation or the gravitational potential on large 
scales within given inflationary models. These, however, are not the physical 
quantities which are observed. One should instead provide a full prediction for 
the second-order radiation transfer function. A preliminary step towards this goal 
has been taken in Ref. [20] (see also [21]) where the full second-order radiation 
transfer function for the CMB anisotropies on large angular scales in a flat uni- 
verse filled with matter and cosmological constant was computed, including the 
second-order generalization of the Sachs-Wolfe effect, both the early and late 
Integrated Sachs-Wolfe (ISW) effects and the contribution of the second-order 
tensor modes. 

There are many sources of NG in CMB anisotropies, beyond the primordial 
one. The most relevant sources are the so-called secondary anisotropies, which 
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arise after the last scattering epoch. These anisotropics can be divided into two 
categories: scattering secondaries, when the CMB photons scatter with electrons 
along the line of sight, and gravitational secondaries when effects are mediated by 
gravity [22]. Among the scattering secondaries we may list the thermal Sunyaev- 
Zeldovich effect, where hot electrons in clusters transfer energy to the CMB 
photons, the kinetic Sunyaev-Zeldovich effect produced by the bulk motion of 
the electrons in clusters, the Ostriker-Vishniac effect, produced by bulk motions 
modulated by linear density perturbations, and effects due to reionization pro- 
cesses. The scattering secondaries are most significant on small angular scales 
as density inhomogeneities, bulk and thermal motions grow and become sizeable 
on small length-scales when structure formation proceeds. 

Gravitational secondaries arise from the change in energy of photons when the 
gravitational potential is time-dependent, the ISW effect, and gravitational lens- 
ing. At late times, when the Universe becomes dominated by the dark energy, 
the gravitational potential on linear scales starts to decay, causing the ISW effect 
mainly on large angular scales. Other secondaries that result from a time depen- 
dent potential are the Rees-Sciama effect, produced during the matter-dominated 
epoch by the time evolution of the potential on non-linear scales. 

The fact that the potential never grows appreciably means that most second 
order effects created by gravitational secondaries are generically small compared 
to those created by scattering ones. However, when a photon propagates from the 
last scattering to us, its path may be deflected because of the gravitational lensing. 
This effect does not create anisotropies, but only modifies existing ones. Since 
photons with large wavenumbers k are lensed over many regions (<~ k /H, where 
H is the Hubble rate) along the line of sight, the corresponding second-order 
effect may be sizeable. The three-point function arising from the correlation of 
the gravitational lensing and ISW effects generated by the matter distribution 
along the line of sight [23, 24] and the Sunyaev-Zeldovich effect [25] are large 
and detectable by Planck [26]. 

Another relevant source of NG comes from the physics operating at the re- 
combination. A naive estimate would tell that these non-linearities are tiny being 
suppressed by an extra power of the gravitational potential. However, the dy- 
namics at recombination is quite involved because all the non-linearities in the 
evolution of the baryon-photon fluid at recombination and the ones coming from 
general relativity should be accounted for. This complicated dynamics might lead 
to unexpected suppressions or enhancements of the NG at recombination. A step 
towards the evaluation of the three-point correlation function has been taken in 
Ref. [27] where some effects were taken into account in the in so-called squeezed 
triangle limit, corresponding to the case when one wavenumber is much smaller 
than the other two and was outside the horizon at recombination. 

These notes, which are based on Refs. [1,2], present the computation of 
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the full system of Boltzmann equations, describing the evolution of the pho- 
ton, baryon and Cold Dark Matter (CDM) fluids, at second order and neglect- 
ing polarization, These equations allow to follow the time evolution of the CMB 
anisotropies at second order on all angular scales from the early epochs, when 
the cosmological perturbations were generated, to the present time, through the 
recombination era. These calculations set the stage for the computation of the 
full second-order radiation transfer function at all scales and for a a generic 
set of initial conditions specifying the level of primordial non-Gaussianity. Of 
course on small angular scales, fully non-linear calculations of specific effects 
like Sunyaev-Zel'dovich, gravitational lensing, etc. would provide a more accu- 
rate estimate of the resulting CMB anisotropy, however, as long as the leading 
contribution to second-order statistics like the bispectrum is concerned, second- 
order perturbation theory suffices. 

These notes are organized as follows. In Section 3 we provide the second- 
order metric and corresponding Einstein equations. In Section 4 the left-hand- 
side of the Boltzmann equation for the photon distribution function is derived 
at second order. The collision term is computed in Section 5. In Section 6 we 
present the second-order Boltzmann equation for the photon brightness function, 
its formal solution with the method of the integration along the line of sight and 
the corresponding hierarchy equations for the multipole moments. Section 7 con- 
tains the derivation of the Boltzmann equations at second order for baryons and 
CDM. Section 8 contains the approximate solution of the Boltzmann equations 
up to first order. Section 9 contains a brief summary of the results. In Appendix 
A we give the explicit form of Einstein's equations up to second order, while in 
Appendix B we provide the first-order solutions of Einstein's equations in various 
cosmological eras. 

In performing the computation presented in these lecture notes, we have mainly 
followed Ref. [28] (in particular chapter 4), where an excellent derivation of 
the Boltzmann equations for the baryon-photon fluid at first order is given, and 
Refs. [1,2] for their second-order extension. Since the derivation at second order 
is straightforward, but lenghty, the reader might benefit from reading the appro- 
priate sections of Ref. [28]. In the Conclusions (Section 9) we have also provided 
a Table which summarizes the many symbols appearing throughout these notes. 



3. Perturbing gravity 

Before tackling the problem of interest - the computation of the Boltzmann equa- 
tions for the baryon-photon and CDM fluids - we first provide the necessary tools 
to deal with perturbed gravity, giving the expressions for the Einstein tensor 
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perturbed up to second order around a flat Friedmann-Robertson- Walker back- 
ground, and the relevant Einstein equations. In the following we will adopt the 
Poisson gauge which eliminates one scalar degree of freedom from the goi com- 
ponent of the metric and one scalar and two vector degrees of freedom from g^. 
We will use a metric of the form 

ds 2 = a 2 (f]) [-e^drj 2 + 2u l dx l dr 1 + (e~ 2 %- + Xl j)dx'dx j ] , (3.1) 

where a(r/) is the scale factor as a function of the conformal time r], and uji 
and Xij are the vector and tensor peturbation modes respectively. Each metric 
perturbation can be expanded into a linear (first-order) and a second-order part, 
as for example, the gravitational potential $ = $W -f $( 2 )/2. However in 
the metric (13. Il l the choice of the exponentials greatly helps in computing the 
relevant expressions, and thus we will always keep them where it is convenient. 
From Eq. ( 13.1b one recovers at linear order the well-known longitudinal gauge 
while at second order, one finds $ (2) = (2) - 2(0 (1) ) 2 and \l/ (2) = ■0 (2) + 
2(?/>W) 2 where (j>^\ ^ and c/> (2 \ V> (2) (with <f>W = $W and = are 
the first and second-order gravitational potentials in the longitudinal (Poisson) 
gauge adopted in Refs. [6, 29] as far as scalar perturbations are concerned. For 
the vector and tensor perturbations, we will neglect linear vector modes since 
they are not produced in standard mechanisms for the generation of cosmological 
perturbations (as inflation), and we also neglect tensor modes at linear order, 
since they give a negligible contribution to second order perturbations. Therefore 
we take cj,; and Xij to be second-order vector and tensor perturbations of the 
metric. 

Let us now give our definitions for the connection coefficients and their ex- 
pressions for the metric (13. U . The number of spatial dimensions is n = 3. Greek 
indices (a, ft, fx, v, ....) run from to 3, while latin indices (a,b, . . . ,i, j, k, . . ., 
m,n, . . .) run from 1 to 3. The space-time metric g^ v has signature (— , +, +, +). 
The connection coefficients are defined as 

r a _ 1 n ap ( d 9 P l , d 9P P dgfr\ 

1 0i - o 9 ^-Tf + -z-z H3T ■ ^-A) 



The Riemann tensor is defined as 

pet pa pet I pet pA pet pA /o o\ 

0pv — 1 0v,p ~ 1 0p.v ' 1 A/i 1 0u ~ 1 Ai/ 1 0n ■ V 3 - 3 ) 

The Ricci tensor is a contraction of the Riemann tensor, R^ v = R a n av and in 
terms of the connection coefficient it is given by 

n fj pet q pa _i_ pa per pa per /o a\ 

*>-p,v — u a 1 fiu U P 1 va ' 1 era 1 pv 1 av 1 pa ■ \ J -^) 
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The Ricci scalar is the trace of the Ricci tensor, R = W 1 .. . The Einstein tensor is 
defined as G M „ = R^ v - ^g^R- 

The Einstein equations are written as G M „ = k 2 T m „, so that k 2 — 87tGn, 
where Gn is the usual Newtonian gravitational constant. 



4. The collisionless Boltzmann equation for photons 

We are interested in the anisotropies in the cosmic distribution of photons and 
inhomogeneities in the matter. Photons are affected by gravity and by Compton 
scattering with free electrons. The latter are tightly coupled to protons. Both 
are, of course, affected by gravity. The metric which determines the gravitational 
forces is influenced by all these components plus CDM (and neutrinos). Our plan 
is to write down Boltzmann equations for the phase-space distributions of each 
species in the Universe. 

The phase-space distribution of particles g(x l , P M ,rj) is a function of spa- 
tial coordinates x i , conformal time r\, and momentum of the particle P M = 
dx* 1 /dX where A parametrizes the particle path. Through the constraint P 2 = 
g^uP^P" = —m 2 , where m is the mass of the particle one can eliminate P° 
and g(x l , P J , 77) gives the number of particles in the differential phase-space vol- 
ume dx 1 dx 2 dx 3 dP 1 dP 2 dP 3 . In the following we will denote the distribution 
function for photons with /. 

The photons' distribution evolves according to the Boltzmann equation 

f = C[/], (4.1) 

where the collision term is due to the scattering of photons off free electrons. 
In the following we will derive the left-hand side of Eq. J4.lt while in the next 
section we will compute the collision term. 

For photons we can impose P 2 = gp, v P li P v — and using the metric (13. U 
in the conformal time r\ we find 



P 2 = a 2 



-e 2 *(P ) 2 + 4+2^P°P < 



= 0, (4.2) 



where we define 

,2 



p z = g VJ P l P3 . (4.3) 
From the constraint (14.21) 



1/2 

P° = e"* f % + 2^P°P 4 ) . (4.4) 
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Notice that we immediately recover the usual zero and first-order relations P° = 
p/a and P° = p(l - $W)/a. 

The components P l are proportional to pn\ where n l is a unit vector with 
n l rii = 8ijn l n3 — 1. We can write P % — Cri 1 , where C is determined by 

<hJ"l" = C 2 a 2 (e- 2 * + X «nV) = p 2 , (4.5) 
so that 

P i = £ n < ( e -a« + Xfcm „*n m )~ 1/2 = 2„< e « f 1 - ± Xfcm „*n m> ) , (4.6) 

where the last equality holds up to second order in the perturbations. Again we 
recover the zero and first-order relations P l = pn l /a and P l = pn l (l + 1 3K 1 >)/a 
respectively. Thus up to second order we can write 

P° = e -»?(1 + Wi n i ) . (4.7) 

Eq. ( 14.6b and (14.7b allow us to replace P° and P l in terms of the variables p and 
n l . Therefore, as it is standard in the literature, from now on we will consider 
the phase-space distribution / as a function of the momentum p = pn l with 
magnitude p and angular direction ri L , f = f(x l ,p, n l , rj). 

Thus, in terms of these variables, the total time derivative of the distribution 
function reads 



(4.8) 



df _ df df dx l df dp df dn 1 
dr) drj dx % drj dp dr\ dn % drj 

In the following we will compute dx % jdr\, dpjdr\ and dn l /drj. 

a) dx % jdr\: 
From 

pi = (W_ = dtfdn = dx^ p0 
dX dr) dX dr\ 

and from Eq. ( |4~6T > and ( R77b 

^ = rte*** (l-uj v? - \ Xkm n k n m ^ . (4. 10) 



b) dp/ dr): 
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For dp/drj we make use of the time component of the geodesic equation 

dP°/d\ = -T° a pP a pP, where d/dX = (drj/dX) d/dq = P° d/dq, and 



dP° 
dq 



-r° 



pap/3 



a/3 pQ > 

Using the metric (13. U we find 



(4.11) 



2T%P a P? = g° 



t dg va _ 5 a/ 3 
dx' 3 dx v 

,2 



a p0 



pap 



2(H + $') (P°) + 4$ iP°P l + 4Ww i P°P i 



2e 



-2* 



/\ -2* j 



1 ./ 



pipj 

(4.12) 



On the other hand the expression (14.71 ) of P° in terms of p and n J gives 

-3— = -~~r e * (l + c^n*) +e * (l+w*n l ) -^-^ 
a?; a dq aq 

p c _^ d{ujin % ) 

a dq 

Thus Eq. ( 14.1 11 1 allows us express dpjd-q as 

1 dp 



p dq 



= -H + - $ i n'e 



/ * 7 



(4.13) 



(4.14) 



where in Eq. ( 14.121 ) we have replaced P° and P % by Eqs. ( 14.71 ) and ( 14.61 ). Notice 
that in order to obtain Eq. ( 14.141 ) we have used the following expressions for the 
total time derivatives of the metric perturbations 



dq 



d<f> d<S> dx i 

drj dx l dq 



k in 



and 



d(o;in l ) 
dq 



duji doji dx 3 
dq dx 3 dr\ 



-n 1 



i9?7 



9a? 



(4.15) 



(4.16) 



where we have taken into account that uj i is already a second-order perturbation 
so that we can neglect dn l /dq which is at least a first order quantity, and we can 
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take the zero-order expression in Eq. ( 14.101 ). dx l /drj = ri. In fact there is also 
an alternative expression for dp /drj which turns out to be useful later and which 
can be obtained by applying once more Eq. ( 14.151 ) 

= -H - ^ +$' + *'- w'i ri - ~X<nV . (4.17) 
p dri dn 2 J 



c) dri 1 /drj: 



We can proceed in a similar way to compute dri 1 /drj. Notice that since in Eq. (14.8b 
it multiplies df /dri 1 which is first order, we need only the first order pertur- 
bation of dri /drj. We use the spatial components of the geodesic equations 

dP l /d\ = -r 4 Q/3 P Q P /3 written as 



dP i 
dr) 



papfi 



a/3 p0 

For the right-hand side we find, up to second order, 



2T l af3 P a P (} = g l 



dxP dx c 



dx v 



= 4(H - y')P l P Q + 2 f 



papP 



pOpk 



(4.18) 



(4.19) 



9$ 



dx 



k r>m 



-2^f-5 km P«P 
ox 1 



«V, + d X \ d Xjk 



dx k dx J 



dxi 



P 3 P 



while the expression ( 14.6b of P % in terms of our variables p and ri in the left-hand 
side of Eq. ( 14.18b brings 



dP i 
drj 



P * 
— e 



dV , a d{p/ a) , 
— ri + — y -^—^-ri 
drj p drj 



dri 
drj 



1 - 2 -Xk m n k n m 



-n e „ 
a 2 



drj 



(4.20) 



Thus, using the expression (14.6b for P % and d4.4b for P° in Eq. ( 14.19b . together 
with the previous result ( 14.141 ). the geodesic equation ( 14.18b gives the following 
expression dri /drj (valid up to first order) 



dri 
drj 



(<f>, k + ^,k)n k ri -*'\ 



(4.21) 
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To proceed further we now expand the distribution function for photons around 
the zero-order value /'°) which is that of a Bose-Einstein distribution 

/ (0) (P, V) = 2 f^—, , (4.22) 

where T(rf) is the average (zero-order) temperature and the factor 2 comes from 
the spin degrees of photons. The perturbed distribution of photons will depend 
also on x 1 and on the propagation direction n l so as to account for inhomo- 
geneities and anisotropies 

f(x i ,p,n i ,r ] ) = fW(p,r)) + f^(x\p,n\ V ) + n\ V ) , (4.23) 

where we split the perturbation of the distribution function into a first and a 
second-order part. The Boltzmann equation up to second order can be written 
in a straightforward way by recalling that the total time derivative of a given i-th 
perturbation, as e.g. df^ /drj is at least a quantity of the i-th order. Thus it is 
easy to realize, looking at Eq. ( 14.8b , that the left-hand side of Boltzmann equation 
can be written up to second order as 

* d f {1} . id f {2) p d 4l±(^ + i^ 



drj drj 2 drj dp drj \ 2 

df (0) d f^m , , T .m , U f21 , 1 



v- 



dp drj 



+ ±$(2) + i^(2) 



- p— — n l - -p— ^r L n l n : > , (4.24) 

op orj 2 op o?7 

where for simplicity in Eq. ( 14.241 ) we have already used the background Boltz- 
mann equation (df / drj) \ ^ = 0. In Eq. (14.24b there are all the terms which will 
give rise to the integrated Sachs-Wolfe effects (corresponding to the terms which 
explicitly depend on the gravitational perturbations), while other effects, such as 
the gravitational lensing, are still hidden in the (second-order part) of the first 
term. In fact in order to obtain Eq. ( 14.24b we just need for the time being to know 
the expression for dp /drj, Eq. ( 14.171 ). 



5. Collision term 

5.1. The Collision Integral 

In this section we focus on the collision term due to Compton scattering 

e(q)7(p) <— > e(q')7(p') , (5.1) 
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where we have indicated the momentum of the photons and electrons involved in 
the collisions. The collision term will be important for small scale anisotropics 
and spectral distortions. The important point to compute the collision term is that 
for the epoch of interest very little energy is transferred. Therefore one can pro- 
ceed by expanding the right hand side of Eq. ( 14. U both in the small perturbation, 
Eq. (14.231 1. and in the small energy transfer. Part of the computation up to second 
order has been done in Refs. [30-32] (see also [33]). In particular Refs. [30, 31] 
are focused on the effects of reionization on the CMB anisotropics thus keeping 
in the collision term those contributions which are relevant for the small-scale ef- 
fects due to reionization and neglecting the effects of the metric perturbations on 
the left-hand side of Eq. ( 14. 11 1. We will mainly follow the formalism of Ref. [31] 
and we will keep all the terms arising from the expansion of the collision term up 
to second order. 

The collision term is given (up to second order) by 

C(p) = C(p)ae* , (5.2) 
where a is the scale factor andQ 

r( s 1 f dq dp' 

° lPJ E(p)J (27r)32£(q) (2 7 r)32 J B(q') (2n)*2E(p<) 

x (2ir) 4 S 4 (q + p-q' -p')\M\ 2 

x { 5 (q')/(p')[l + /(p)]-5(q)/(p)[l + /(p')]} (53) 

where E(q) = (q 2 + m 2 ) 1 / 2 , M is the amplitude of the scattering process, 
S 4 (q + p-q'-p') = £ 3 (q + p-q' - p')S(E(q) +p - E(q') -p') ensures the 
energy-momentum conservation and g is the distribution function for electrons. 
The Pauli suppression factors (1 — g) have been dropped since for the epoch 
of interest the density of electrons n e is low. The electrons are kept in thermal 
equilibrium by Coulomb interactions with protons and they are non-relativistic, 
thus we can take a Maxwell-Boltzmann distribution around some bulk velocity v 

( 2 ^ \ 3/2 f (q-™,v) 2 l 

\m e l e J I 2m e J-e J 

By using the three dimensional delta function the energy transfer is given by 
E(q) — _E(q + p — p') and it turns out to be small compared to the typical 



1 The reason why we write the collision term as in Eq. 15. 2t is that the starting point of the Boltz- 
mann equation requires differentiation with respect to an affine parameter A, df /d\ = C In mov- 
ing to the conformal time rj one rewrites the Boltzmann equation as df /drj = C'(P°)~ , with 
P° = drj/dX given by Eq. j4.7l . Taking into account that the collision term is at least of first order, 
Eq. 15. 2\ then follows. 
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thermal energies 



£(q)-£(q + p-p') 



(P- 



= 0(Tq/m e ) . 



(5.5) 



In Eq. (15.5b we have used E(q) = m e + q 2 /2m e and the fact that, since the 
scattering is almost elastic (p ~ p'), (p — p') is of order p ~ T, with g much 
bigger than (p — p'). In general, the electron momentum has two contributions, 
the bulk velocity (q = m e v) and the thermal motion (q ~ (nifT) 1 / 2 ) and thus 
the parameter expansion q/m e includes the small bulk velocity v and the ratio 
(I/me) 1 / 2 which is small because the electrons are non-relativistic. 

The expansion of all the quantities entering the collision term in the energy 
transfer parameter and the integration over the momenta q and q' is described in 
details in Ref. [31]. It is easy to realize that we just need the scattering amplitude 
up to first order since at zero order .g(q') = g(q + p — p') = g(q) and S(E(q) + 
p — -B(q') —p') = S(p — p') so that all the zero-order quantities multiplying \M\ 
vanish. To first order 



\M\ 2 = 6Tra T m 2 e [(l + cos 2 6>) - 2cos6>(l - cos6>)q • (p + p')/r 



(5.6) 



where cos# — n • n' is the scattering angle and <tt the Thompson cross-section. 
The resulting collision term up to second order is given by [31] 



3n e a T f , ,dfl' 



4 2 >( p ,pO + c£W) + c£W) 



.(1) 



(p,p')+4 2) (p,p') + 4 2) (p,p') 



(2), 



(5.7) 



where we arrange the different contributions following Ref. [31]. The first order 
term reads 



c (1) (P,P') - (1 + cos 2 



S(p- P 'Kf {1) (p')-f (1, (p)) 



(i)/ 



+ (/ (0) (P')-/ (0) (P))(P-P')-V^|^ 



(5.8) 



while the second-order terms have been separated into four parts. There is the 
so-called anisotropy suppression term 

4 2) (P, P') = \ (1 + cos 2 0) S(p - p')(/ (2) (p') - / (2) (P)) ; (5.9) 



Cosmic Microwave Background Anisotropics up to Second Order 



17 



a term which depends on the second-order velocity perturbation defined by the 
expansion of the bulk flow as v = + v^ 2 ' /2 

ci 2 )(p,p') = I(l+Cos^)(/( )( P ')-/ (0) (P))(P-P')-V (2) dS % pl P,) ;(5-10) 
a set of terms coupling the photon perturbation to the velocity 

cS(P-P') = (/ (1) (P')-/ (1) (P)) (l + cos 2 0)(p-p')-v 

x 96{j? d ~ ^ - 2cosfl(l - cos6)S(p - p')(n + n') ■ v 



and a set of source terms quadratic in the velocity 

c&W) = (/ (0) (P')-/ (0) (P)) (P-P')-V 
(p-pi)-vd 2 5(p-p>) 



(1 + cos 2 6») 



2 <9p' 2 
2cos6»(l -cos#)(n + n') v 



d6 (p-p') 
dp' 



(5.11) 



The last contribution are the Kompaneets terms describing spectral distortions to 
theCMB 



4\p,p') = (l + cos 2 0) 

d 2 S(p-p') 
x dp 12 



2M (P-P') 



2m 



+ 



(/( )(p')-/(°)(p))T e (5.12) 

(/ (0) (p') + / (0) (p) + 2.f (0) (p')/ (0) (rt) 

2(p - f>')cos0(l - cos 2 6») 



m e 



x / (0) b')(l + / (0) (p))(/ (0) b')-/ (0) (p) 



*(p-j/) 

9(5 (p — p') 
dp~' 



Let us make a couple of comments about the various contributions to the col- 

(2) 

lision term. First, notice the term c v (p, p ) due to second-order perturbations 
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in the velocity of electrons which is absent in Ref. [31]. In standard cosmo- 
logical scenarios (like inflation) vector perturbations are not generated at linear 
order, so that linear velocities are irrotational i/ 1 '* = d i v^ 1 '. However at sec- 
ond order vector perturbations are generated after horizon crossing as non-linear 
combinations of primordial scalar modes. Thus we must take into account also a 
transverse (divergence-free) component, i/ 2 ) 1 = d l v^ + v£ with diVrp * = 0. 
As we will see such vector perturbations will break azimuthal symmetry of the 
collision term with respect to a given mode k, which instead usually holds at 
linear order. Secondly, notice that the number density of electrons appearing in 
Eq. ( 15.7b must be expanded as n e = n e {\ + 5 e ) and then 

4 X) C (1) (P,P') (5.13) 

gives rise to second-order contributions in addition to the list above, where we 
split S e — 5e + 8^/2 into a first- and second-order part. In particular the 
combination with the term proportional to v in c^(p, p') gives rise to the so- 
called Vishniac effect, as discussed in Ref. [31], 

5.2. Computation of different contributions to the collision term 

In the integral ( 15. 7t over the momentum p' the first-order term gives the usual 
collision term 



C^(p) = n e a T 



/i 1) ( P ) + ^(v.n)-/«-^v.n 



(5.14) 



where one uses the decomposition in Legendre polynomials 

/( 1 )(x,p,n)=^(2^ + l)/i 1) (p)P,(cos 1 ?), (5.15) 
i 

where -d is the polar angle of n, cosi? = n • v. 

In the following we compute the second-order collision term separately for 
the different contributions, using the notation C(p) = C' 1 ' (p) + (p)/2. We 
have not reported the details of the calculation of the first-order term because for 

(2) (2) 

its second-order analog, c A (p, p') + c v (p, p'), the procedure is the same. The 
important difference is that the second-order velocity term includes a vector part, 
and this leads to a generic angular decomposition of the distribution function (for 
simplicity drop the time dependence) 



I m=-l " 
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such that 



/£ = (-<)- 



21 + 1 

47T 



dn/«y/ ro (n). 



(5.17) 



Such a decomposition holds also in Fourier space [34]. The notation at this stage 
is a bit confusing, so let us restate it: superscripts denote the order of the per- 
turbation; the subscripts refer to the moments of the distribution. Indeed at first 
order one can drop the dependence on m setting m = using the fact that the 
distribution function does not depend on the azimuthal angle <f>. In this case the 
relation with fP~' is 

a)4 2) (p, P '): 

The integral over p' yields 



(5.18) 



4p 
dQ' 

ill 



[l + (n.n') 2 ][/ (2) (p')-/ (2, (p)] 



dp'p'6(p — p') 

(5.19) 



To perform the angular integral we write the angular dependence on the scattering 
angle cos6* = n • n' in terms of the Legendre polynomials 



[l + (n.n') 2 ] = - 



l + ip 2 (n.n') 



i + jE ^ 2m (n)r 2 ;„(n') 



4tt 



m=-2 



21+ 1 



(5.20) 



where in the last step we used the addition theorem for spherical harmonics 

2 



Pi = 



47T 



(5.21) 



Using the decomposition ( 15.17b and the orthonormality of the spherical harmon- 
ics we find 



Ca ) (P) = n e a T 



(5.22) 
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It is easy to recover the result for the corresponding first-order contribution in 
Eq. ( ETBb by using Eq. ( fgTTgb . 

b)ci 2) (p,p'): 

Let us fix for simplicity our coordinates such that the polar angle of n' is 
defined by p! = v( 2 ) • n' with <f>' the corresponding azimuthal angle. The contri- 

(2) 

bution of c v (p, p') to the collision term is then 

x |* - P ' f /) fj g[l + (p • p') 2 ] ■ (5.23) 

We can use Eq. (15.20b which in our coordinate system reads 



2 ^ (2 + m)\ 2 K ' 2 y ' 

m— — 2 K ' 



(5.24) 



4 
3 

so that 

J ^P 2 (n • n') = P 2 (n • V ( 2 >)P 2 (n' • v< 2 >) = P 2 (fx)P 2 (»') . (5.25) 

By using the orthonormality of the Legendre polynomials and integrating by parts 
over p' we find 

Ci 2) (p) = -n e a T p^v( 2 ) ■ n . (5.26) 
op 

As it is clear by the presence of the scalar product v^ 2 ) • p the final result is 
independent of the coordinates chosen. 

c)4 2 i(p,p'): 

Let us consider the contribution from the first term 



4Vp=p')= (i + cos 2 fl) (/«(p')-/ (1) (p)) (p-p')-v d5{p ~ p,) 



dp' 

where the velocity has to be considered at first order. In the integral ( 15. 7t it brings 



1^(2) _ 3n e a T 



2 At, « 4p 

x - pV) 



"-Jdp'p'^^J^Hp')-^)] 

p2-K Jjj 

J ^(l + cos 2 0), (5.27) 
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The procedure to do the integral is the same as above. We use the same relations 
as in Eqs. ( 15.24b and ( 15.25b where now the angles are those taken with respect to 
the first-order velocity. This eliminates the integral over cf>', and integrating by 
parts over p' yields 



3n e o~TV f 1 dfi' 



(5.28) 



<9p 

We now use the decomposition ( 15.15b and the orthonormality of the Legendre 
polynomials to find 

/ ^m7 (1 W)W) = E / ^mW)W)/J x) (p) 



^P l(M ') + ^P3(/x') 



= |/X (1) (P) + ^3 (1) (P) 



(5.29) 



where we have used fi'P 2 (fi')P l (fjt') = §Pi(//) + tPsO'), with P x (//) = //. 



Thus from Eq. (15.28b we get 



/ (1) (p)-/o (1) (p)- 



dp 



-P 2 (^-n)[f^\p)+p 



,d/ 2 (1) (^ 
dp 



P 



dp 



dp 



(5.30) 



In c*- 2 )(p, p') there is a second term 

c (2) 

whose contribution to the collision term is 

2>n e (jTV 



2 C 'a](/7)(p) 



= -2cos0(l - cos<9) (/ (1) (p') - / (1) (P)) S(p-p')(n + n') • v, 
J dp'p'5{p-p') £^V)(p') 



2p 
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f W (j>))(f* + n')J ^ CO s#(l - costf) . (5.31) 



This integration proceeds through the same steps as for C^K^ (p). In particular 
by noting that cos6»(l - cos6>) = -1/3 + Pi(cos0) - 2P 3 (cos6»)/3, Eqs. (15341 
and (I5.25l l allows to compute 



J ^cosfl(i - cose) = -1 + f\o*)PiOO - , 

and using the decomposition ( 15.15b we arrive at 



(5.32) 



1 



C 



(2) 



2 Av(II) 



(P) 



+ u 



, e a T |v-n/ 2 (1) (p)(l-P 2 (v-n)) 
ip 2 (v-n) (3A (1) (p)-3/3 (1) (p) 



(5.33) 



We then obtain 



~^(p) 



= n e er T <v-n 



/ (1) (p)-/o (1) (P)-P^^-/ 2 (1) W 



g/ 2 (1) (p) 
dp 



~P 2 (v-n)(/ 2 (1 V) -j, 

a/i (1) (p) , c ,(d, . , 3 a/^fr) 



2A (1) (p) 



(5.34) 



As far as the remaining terms, these have already been computed in Ref. [31] 
(see also Ref. [30]) and here we just report them 

d)cS?(p,p'): 

The term proportional to the velocity squared yield a contribution to the colli- 
sion term 



1 



C#(P) = 



"eOT j(v • n) 2 

d f(0) 



ii 2 3 2 / (0) 



dp 2(T dp 2 
2 5 2 /W 



<9p 2CT dp 2 



(5.35) 



e) c^(p,p'): 
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~C%\p) 



1 d J 4 
= — — <p 



m\ dp \ 



T e ^ + /(°)(l + /(°)) 
op 
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(5.36) 



Finally, we write also the part of the collision term coming from Eq. d5.13t 



-/ 2 (1) P 2 (v.n)-/ 



(i) 



/ (1) (P) 
dp 



(5.37) 



5.5. Fj/w/ expression for the collision term 

Summing all the terms we find the final expression for the collision term (15.7b up 
to second order 



C(p) = C«(p) + ic( 2 >(p) 



with 

C«(p) = n e a T 
and 



/o ( %) + iU (1) iM*-n)-/«- 



dp 



-v • n 



(5.38) 



(5.39) 



^ 2 )( P ) = n^i/W&O-i £ ^/SWMn) 

771— — 2 



5/ (2) (p)+4 1) [/o (1) (F) + ^/ 2 (1) ^(v-n)-/W 



p— — v • n 

op 



n + v • n 



/ (1) (P) 



x I/ 2 (1) (P)-P 



dp 
dp 

iP 2 (v-n)( -/! (1) (p)+p 



/r>)+2^-n) 



2A (1) (p)+p 



,fl /i (1) (p) 

dp 



3 A/j"(p) 
2^ dp 



+ (v • n) 2 



ii 2 a 2 /(°) 



dp ' 20 P dp 2 
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P- 



dp 



J_d_ 

m 2 dp 



P 



20 P 



dp 2 



(5.40) 



Notice that there is an internal hierarchy, with terms which do not depend on the 
baryon velocity v, terms proportional to v ■ n and then to (v • n) 2 , v and v 2 
(apart from the Kompaneets terms). In particular notice the term proportional to 
S^v ■ n is the one corresponding to the Vishniac effect. We point out that we 
have kept all the terms up to second order in the collision term. In Refs. [30, 31] 

(2) 

many terms coming from c Av have been dropped mainly because these terms are 
proportional to the photon distribution function / W which on very small scales 
(those of interest for reionization) is suppressed by the diffusion damping. Here 
we want to be completely general and we have to keep them. 



6. The Brightness equation 



6.1. First order 



The Boltzmann equation for photons is obtained by combining Eq. ( 14.24b with 
Eqs. d5.39t - d5.40b . At first order the left-hand side reads 



dq 



dfW 0/(0) 0$(1) dx i 0/(°) 0*W 



dp dx % dt] 



dp drj 



(6.1) 



At first order it is useful to characterize the perturbations to the Bose-Einstein 
distribution function d4.22t in terms of a perturbation to the temperature as 



f(x\p,n\ V ) = 2 
Thus it turns out that 



exp 



T(r?)(l + e«) 



i -i 



- 1 



(i) 



dp 



(6.2) 



(6.3) 



where we have used the fact that <9//90|e=o = —pdf^°'/dp. In terms of this 
variable 0^ the linear collision term d5.39l l will now become proportional to 
—pdf(°> I dp which contains the only explicit dependence on p, and the same 
happens for the left-hand side, Eq. d6.U . This is telling us that at first order GW 
does not depend on p but only on x l , n l , rj, = (x\ n l , r). This is well 
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known and the physical reason is that at linear order there is no energy transfer 
in Compton collisions between photons and electrons. Therefore, the Boltzmann 
equation for 0W reads 



dew 

= n e o~T a 



dx l 



drj 



9^ J P 2 (v-n) -9 ( 



(6.4) 



where we made us of /j 1 ' = —pdf^/dpQg, according to the decomposition 
of Eq. (15.15b . and we have taken the zero-order expressions for dx % /drj, dropping 
the contribution from dn l jdr\ in Eq. ( 14.81 ) since it is already first-order. 

Notice that, since is independent of p, it is equivalent to consider the 
quantity 



A«(a 



Jdpp 3 fM 
J dpp 3 fW ' 



(6.5) 



being A' 1 ) = 40W at this order. The physical meaning of A' 1 ) is that of a 
fractional energy perturbation (in a given direction). From Eq. (14.241) another 
way to write an equation for A^ - the so-called brightness equation - is 



d_ 

drj 



(1) JL 



n e o~T a 



_ 4 £ 
or] \ 

A< 1) P 2 (v-n) 



1 



4v ■ n 



(6.6) 



6.2. Second order 

The previous results show that at linear order the photon distribution function 
has a Planck spectrum with the temperature that at any point depends on the 
photon direction. At second order one could characterize the perturbed photon 
distribution function in a similar way as in Eq. (16.21 ) 



f(x\p,n\r]) = 2 



exp 



T(r)) e c " 



- 1 



(6.7) 



where by expanding = 0W + 0( 2 )/2 + ... as usual one recovers the first- 
order expression. For example, in terms of 0, the perturbation of is given 
by Eq. ( 16.31 ). while at second order 



/(2) _ pdf (0) 1 

2 2 dp 2 



,a 2 /(°) 

dp 2 



(6.8) 
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However, as discussed in details in Refs. [30, 31], now the second-order per- 
turbation 0( 2 ) will not be momentum independent because the collision term 
in the equation for O' 2 ' does depend explicitly on p (defining the combination 
— (pS/^/dp) -1 /' 2 ' does not lead to a second-order momentum independent 
equation as above). Such dependence is evident, for example, in the terms of 
C^(p), Eq. (15.40b . proportional to v or v 2 , and in the Kompaneets terms. The 
physical reason is that at the non-linear level photons and electrons do exchange 
energy during Compton collisions. As a consequence spectral distorsions are 
generated. For example, in the isotropic limit, only the Kompaneets terms sur- 
vive giving rise to the Sunyaev-Zeldovich distorsions. As discussed in Ref. [30], 
the Sunyaev-Zeldovich distorsions can also be obtained with the correct coeffi- 
cients by replacing the average over the direction electron (v 2 ) with the mean 
squared thermal velocity (v 2 h ) — 3T e /m e in Eq. ( 15.401 ). This is due simply to 
the fact that the distinction between thermal and bulk velocity of the electrons 
is just for convenience. This fact also shows that spectral distorsions due to the 
bulk flow (kinetic Sunyaev-Zeldovich) has the same form as the thermal effect. 
Thus spectral distorsions can be in general described by a global Compton y- 
parameter (see Ref. [30] for a full discussion of spectral distorsions). However 
in the following we will not be interested in the frequency dependence but only 
in the anisotropies of the radiation distribution. Therefore we can integrate over 
the momentum p and define [30, 3 1 ] 

5 ' J dpp3 /(0) ' { ' 

as in Eq. ( 16.5b . 

Integration overp of Eqs. rt4.24b -( f5~40T > is straightforward using the following 
relations 

dpp 3 p^-— = -AN; / dpp 3 p 2 -^-5- = 20 7V; 



dp J dp 2 



Jdpp 3 f w = NA (1) ; J dpp 3 p^- = -4NA (1) . (6.10) 

Here N — J dpp 3 f^ is the normalization factor (it is just proportional the 
background energy density of photons p 7 ). At first order one recovers Eq. ( 16.6l l. 
At second order we find 



1 d 
2~ch] 



A( 2 ) +4$< 2 ) 



d 

df] 

2 A f m + & 2 A + 4^V + 2%^nV 
or] V J orj orj 



4A« f*(D'_$W, 
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.(1) 



l-A i 2 1) P 2 (v • n) - A (1) + 4v • n ) + 4v (2) ■ n 



+2(v • n) 



vA^ (4 + 2P 2 (v • n)) + 14(v • n) z - 2v 2 



A«+3A«-A« l--P 2 (v-n 



where we have expanded the angular dependence of A as in Eq. ( 15.161 1 



(6.11) 



A«(x,n)=£ £ A ^(x)H)' 



47T 

2^ + 1 



^m(n) , 



(6.12) 



with 



H)" 



^y^/dOA^n), 



(6.13) 



where we recall that the superscript stands by the order of the perturbation. At 
first order one can drop the dependence on m setting m = so that A^ = 
(— i)~ e (2£+ 1)6 mo A^' . In Eq. ( 16.111 ) we have introduced the differential optical 
depth 



t = ~n e o-T& ■ 



(6.14) 



It is understood that on the left-hand side of Eq. ( 16.111 ) one has to pick up for the 
total time derivatives only those terms which contribute to second order. Thus 
we have to take 



ld_ 
2~ch] 



A< 2 > +4$( 2 



d 
dr) 



A«+4$« 



(2) 



(6.15) 



2 \dr] dx l J V / v ; 

L » » J an 1 

where we used Eqs. ( 14.101 ) and ( 14.211 ). Notice that we can write dA^ 1 ' /dn l — 
(SAW /dx i )(dx l /dri l ) = {dA^> / dx l ){ri — rji), from the integration in time of 
Eq. ( 14.101 ) at zero-order when n l is constant in time. 
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6.3. Hierarchy equations for multipole moments 

Let us now move to Fourier space. In the following, for a given wave-vector k 
we will choose the coordinate system such that e3 = k and the polar angle of the 
photon momentum is d, with /i = cos?9 = k ■ n. Then Eq. (16. lit can be written 

as 



A< 2 >' + ik^ _ T ' A (2) = S ( k) n ^ ( 



(6.16) 



where S(k, n, 77) can be easily read off Eq. (16.1 lb . We now expand the temper- 

(2) 

ature anisotropy in the multipole moments A}' in order to obtain a system of 
coupled differential equations. By applying the angular integral of Eq. ( 16.131 ) to 
Eq. ( 16.161 ) we find 



A£>'(M) = fc 



K lm . (2) _ Kl+l,m . (2) 

2i _i e-i, m 2 e + 3 £+1 ' m 



T'A?2 + S tm (6.17) 



where the expansion coefficients of the source term are given by 



S, 



j±i<>ei 



2t' 



d z k x 



^(^^(kajka ■ k! + (*W(k x ) + $ (1) (ki)) 



x A«(k 2 )-i^(k 1 )AW(k 



16fe 



d 3 /ci 
(2^)3 



$(i)(k 1 )$( 1 )(k 



8r'[(<5« + $)v] m + 6r'(A«v) ro _ 2 r'(A 2 1) v 



(xIf(l)V$( 1 )) r 



£<i + t' 



Ski 

X 5i2 S m + 



(5«(kx) + $( 1 )(k 1 ))A«(k 2 ) - *- v (k 1 )Ai 1 ) (k 2 ) 



(2tt) 3 

(??-%)(* (1) + < f >(1) )(ki)k 1 -k 2 
(-1)— (2* + l)£ £ (2^ + 1) 



8*«(k!) + 2r'( ( 5«(k 1 ) + $M(k!)) 



<" m'— - 



8A^V$ (1) - 2($W 
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t'zH)^(-1)- 1 (2£+1)£ ^ (2f'+l) 

m' = -l 

f ' 1 £ \ / £" 1 £ 
) \ m' -m 

14r'H)^(-l)- ro £ ' 1 



5^ 2 A^v 



47T. 



-Y 



' 1 * m /(ti)(*ir m «(k)-fciir m «(ti) 



< ; (ki) 



fc 2 



1 1 £ 




2^+1 

47T 



m' m" -m )+ 2 fo-*)H)" 
l 



£ E 



L m' ,m" 



d 3 h 



J (2tt) 3 V 2L+ 1 V 3 



47T / 47T 



x (-1)" 

x A«(k x )($ (1) + * (1) )(k 2 )fc 1 r 1 * m ,(k 1 ) (fcy; m „(k) 



fcxYi^^kO / dfiy lm , (n)yi m »(n)y L0 (n)y,_ m (n), (6.18) 



where k 2 = k — ki and fc 2 = |k 2 |. In Eq. ( 16.181 ) it is understood that \m\ < I. 

Let us explain the notations we have adopted in writing Eq. ( 16.181 ). The baryon 
velocity at linear order is irrotational, meaning that it is the gradient of a potential, 
and thus in Fourier space it is parallel to k, and following the same notation of 
Ref. [37], we write 



vW(k) = -z^ iJ (k)k 



(6.19) 



The second-order velocity perturbation will contain a transverse (divergence- 
free) part whose components are orthogonal to k = e3, and we can write 



r( 2 )(k) = -w<j 2) (k)e 3 + u 



(2) 



m=±l 



e 2 T »ei 
V2 



(6.20) 



where ej form an orhtonormal basis with k. The second-order perturbation uii 
is decomposed in a similar way, with u>±\ the corresponding components (in 
this case in the Poisson gauge there is no scalar component). Similarly for the 
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tensor perturbation Xij we have indicated its amplitudes as x±2 in the decompo- 
sition [36] 



Xij 



^2 -YgXm(ei±ie 2 )i(ei±ie 2 )j 



(6.21) 



m=±2 



We have taken into account that in the gravitational part of the Boltzmann equa- 
tion and in the collsion term there are some terms, like <s£ v> which still can be 
decomposed in the scalar and transverse parts in Fourier space as in Eq. ( 16.201 1. 
For a generic quantity /(x)v we have indicated the corresponding scalar and 
vortical components with (/v) m and their explicit expression is easily found by 
projecting the Fouriermodes of /(x)v along the k = e3 and (e2=R e i) directions 



(/v) m (k) 



d 3 h 
(2tt) 3 '' 



ki)/(k 2 )y 1 * m (k 1 ) A /— . 



Similarly for a term like /(x)Vg(x) we used the notation 

(/V<?) m (k) = -J ^fc lff (k 1 )/(k 2 )Y 1 * ro (k 1 )^ 



(6.22) 



(6.23) 



Finally, the first term on the right-hand side of Eq. ( 16.17) has been obtained by 
using the relation 



zk-nA^k) = ^A^(k) 



Cm 



im 



21+1 
K lm ^(2) 



G lm , (6.24) 



Kim . (2) 

21- \ -t- l - m ' 2£ + 3 e+1 ' m 



where Gg m = (— i) 1 \J j (21 + l)l^ m (n) is the angular mode for the decom- 
position ( 16.12) and Kim = VP — vp?. This relation has been discussed in 
Refs. [34, 36] and corresponds to the term n 4 <9A( 2 ) / 'dx l in Eq. ( 16.1 1) . 

As expected, at second order we recover some intrinsic effects which are 
characteristic of the linear regime. In Eq. (16.17) the relation ( 16.24) represents 
the free streaming effect: when the radiation undergoes free-streaming, the in- 
homogeneities of the photon distribution are seen by the obsever as angular 
anisotropics. At first order it is responsible for the hierarchy of Boltzmann 
equations coupling the different £ modes, and it represents a projection effect 
of fluctuations on a scale k onto the angular scale i ~ kr\. The term rAl' 
causes an exponential suppression of anisotropies in the absence of the source 
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term Sg m . The first line of the source term ( 16.18b just reproduces the expres- 
sion of the first order case. Of course the dynamics of the second-order metric 
and baryon-velocity perturbations which appear will be different and governed 
by the second-order Einstein equations and continuity equations. The remaining 
terms in the source are second-order effects generated as non-linear combinations 
of the primordial (first-order) perturbations. We have ordered them according to 
the increasing number of I modes they contribute to. Notice in particular that 
they involve the first-order anisotropies A^ and as a consequence such terms 
contribute to generate the hierarchy of equations (apart from the free-streaming 
effect). The source term contains additional scattering processes and gravita- 
tional effects. On large scales (above the horizon at recombination) we can say 
that the main effects are due to gravity, and they include the Sachs-Wolfe and the 
(late and early) Sachs-Wolfe effect due to the redshift photons suffer when trav- 
elling through the second-order gravitational potentials. These, toghether with 
the contribution due to the second-order tensor modes, have been already studied 
in details in Ref . [20] . Another important gravitational effect is that of lensing 
of photons as they travel from the last scattering surface to us. A contribution of 
this type is given by the last term of Eq. ( 16.181 ). 

6.4. Integral solution of the second-order Boltzmann equation 

As in linear theory, one can derive an integral solution of the Boltzmann equa- 
tion ( 16.1 U in terms of the source term S. Following the standard procedure 
(see e.g. Ref. [28, 35]) for linear perturbations, we write the left-hand side as 
A< 2 )' + ikfiAW - t'A^ = e - lkm+T d[A^e lk ^- T }/dri in order to derive the 
integral solution 



where rjo stands by the present time. The expression of the photon moments A^i 
can be obtained as usual from Eq. ( 16.131 ). In the previous section we have already 
found the coefficients for the decomposition of source term S 



In Eq. ( 16.251 ) there is an additional angular dependence in the exponential. It is 
easy to take it into account by recalling that 




(6.25) 




(6.26) 




(6.27) 
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Thus the angular integral ( 16.13b is computed by using the decomposition of the 
source term ( 16.26b and Eq. ( 16.27b 

fVo 

A^(Mo) = (-l)- m H)- £ (2^+l) / d V e- T W 

Jo 

1% 

X Yj (- l )' 2S '2m ! Xi i ' 1 ^( fc ( , )-')o)) 

£ 2 m 2 = -e 2 h 

where the Wigner 3 — j symbols appear because of the Gaunt integrals 

C~ = J d 2 m imi (n)Y hm2 (h)Y hm3 (n) 
(2h + 1) {2l 2 + 1) (2l 3 +iy 

47T 

(6.29) 



I V mi m2 TO3 

Since the second of the Wigner 3-j symbols in Eq. ( 16.28b is nonzero only if 
m = ni2, our solution can be rewritten to recover the corresponding expression 
found for linear anisotropies in Refs. [34, 36] 

Ai §t+f - fj J2 Se 2mJ t m) lHm ~ V)} , (6.30) 

where jf 2,m \k{rjQ — rf)] are the so called radial functions. Of course the main 
information at second order is included in the source term containing different ef- 
fects due to the non-linearity of the perturbations. In the total angular momentum 
method of Refs. [34, 36] Eq. ( 16.30b is interpreted just as the intergration over the 
radial coordinate (x = rjo — if) of the projected source term. Another important 
comment is that, as in linear theory, the integral solution ( 16.28b is in fact just a 
formal solution, since the source term S contains itself the second-order photon 
moments up to / = 2 (see Eq. ( 16.18b ). This means that one has anyway to resort 
to the hierarchy equations for photons, Eq. ( 16.17b , to solve for these moments. 
Neverthless, as in linear theory [35], one expects to need just a few moments 
beyond £ = 2 in the hierarchy equations, and once the moments entering in the 
source function are computed the higher moments are obtained from the integral 
solution. Thus the integral solution should in fact be more advantageous than 
solving the system of coupled equations ( 16.17b . 
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7. The Boltzmann equation for baryons and cold dark matter 

In this section we will derive the Boltzmann equation for massive particles, which 
is the case of interest for baryons and dark matter. These equations are necessary 
to find the time evolution of number densities and velocities of the baryon fluid 
which appear in the brightness equation, thus allowing to close the system of 
equations. Let us start from the baryon component. Electrons are tightly coupled 
to protons via Coulomb interactions. This forces the relative energy density con- 
trasts and the velocities to a common value, S e = S p = Sb and v e = v p = v, so 
that we can identify electrons and protons collectively as "baryonic" matter. 

To derive the Boltzmann equation for baryons let us first focus on the colli- 
sionless equation and compute therefore dg/drj, where g is the distribution func- 
tion for a massive species with mass m. One of the differences with respect to 
photons is just that baryons are non-relativistic for the epochs of interest. Thus 
the first step is to generalize the formulae in Section 4 up to Eq. ( 14.211 ) to the case 
of a massive particle. In this case one enforces the constraint Q 2 = g^uQ^Q" = 
—m 2 and it also useful to use the particle energy E = y/ q 2 + to 2 , where q is 
defined as in Eq. (I4.3t . Moreover in this case it is very convenient to take the 
distribution function as a function of the variables q l = qn l , the position x l and 
time ?/, without using the explicit splitting into the magnitude of the momentum 
q (or the energy E) and its direction n % . Thus the total time derivative of the 
distribution functions reads 

dq dq dq dx l dq dq l 

— = — H — 1 -— . (7.1) 

drj drj dx l dr\ dq 1 dr\ 

We will not give the details of the calculation since we just need to replicate the 
same computation we did for the photons. For the four-momentum of the particle 
notice that Q l has the same form as Eq. ( 14.61 1, while for Q° we find 



. a E[l + ^-y (7.2) 

In the following we give the expressions for dx l /drj and dq 1 /drj, 

a) As in Eq. ( 14.101 ) dx l /dr) = Q l /Q° and it turns out to be 

f = |»'e«(l-^f)(l4»,„»V»). ,7.3, 

b) For dq 1 /d-q we need the expression of Q % which is the same as that of Eq. ( 14.61 ) 
Q i = -e* ( 1 - lxk m n k n m ) . (7.4) 



a \ 2' 
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The spatial component of the geodesic equation, up to second order, reads 



dQ 



1 



-2- = -2(W-*0 In'e 



+ <4 -<<)] + 



WwMjfe - -(x*,-,fc + x\j - X 3 k) 



q j q k 
Ea 



(7.5) 



Proceeding as in the massless case we now take the total time derivative of 
Eq. ( 17. 4t and using Eq. ( 17.51 ) we find 

drj ' E E 



E(^' + H^)-{x\+^ k -ui)Eq k 



+ 



q j q k 
~E~ 



(7.6) 



We can now write the total time derivative of the distribution function as 

1 



dg 
dr) 



or} E 



1 — UJiU 
k 



-Xk m n k n m 



-{H - *V + * fc ^-e*+* - ^Ee*^ - * Ae*+* 
E E 



EW + H^)-{x l ' k + u>l k -^)Eq 
1 
2 



k 



+ { Hio l 5 3k - -(x\. k + X\j ~ Xjl) i -g 



q 3 q k 



dg_ 

dq l 



(7.7) 



This equation is completely general since we have just solved for the kinematics 
of massive particles. As far as the collision terms are concerned, for the system 
of electrons and protons we consider the Coulomb scattering processes between 
the electrons and protons and the Compton scatterings between photons and elec- 
trons 



^(x.q,„) 

f(x,Q„) 



( c ep)QQ'q' ~\~ (Ce~{)pp'q' 
( c ep)qq'Q' i 



(7.8) 
(7.9) 



where we have adopted the same formalism of Ref. [28] with p and p' the initial 
and final momenta of the photons, q and q' the corresponding quantities for the 
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electrons and for protons Q and Q'. The integral over different momenta is 
indicated by 

d 3 p f d 3 p' f d 3 q' 

' 1 1 1 , (7.10) 



(2tt) 3 J (2tt) 3 J (2tt) 3 

and thus one can read c ej as the unintegrated part of Eq. d5.2t , and similarly for 
c ep (with the appropriate amplitude |M| 2 ). In Eq. ( 17.8b Compton scatterings be- 
tween protons and photons can be safely neglected because the amplitude of this 
process has a much smaller amplitude than Compton scatterings with electrons 
being weighted by the inverse squared mass of the particles. 

At this point for the photons we considered the perturbations around the zero- 
order Bose-Einstein distribution function (which are the unknown quantities). 
For the electrons (and protons) we can take the thermal distribution described by 
Eq. ( 15.41 ). Moreover we will take the moments of Eqs. ( I7.8l )-( |7!9] l in order to find 
the energy-momentum continuity equations. 

7.1. Energy continuity equations 

We now integrate Eq. ( 17.71) over d 3 q/(2ir) 3 . Let us recall that in terms of the 
distribution function the number density n e and the bulk velocity v are given by 

(7.11) 



(27T) 



and 

1 f d 3 q qn l 

V =-Jwr g ^' (7 - 12) 

where one can set E ~ m e since we are considering non-relativistic particles. 
We will also make use of the following relations when integrating over the solid 
angle dfi 

J dfln 1 = J dQ n l n ] n k = , J nV = . (7.13) 

Finally notice that dE/dq = q/E and dg/dq = (q/E)dg/dE. 

Thus the first two integrals just brings n' e and (ra e v'),;. Notice that all the 
terms proportional to the second-order vector and tensor perturbations of the met- 
ric give a vanishing contribution at second order since in this case we can take 
the zero-order distribution functions which depends only on 77 and E, integrate 
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over the direction and use the fact that S^Xij — 0- The trick to solve the remain- 
ing integrals is an integration by parts over <?. We have an integral like (the one 
multiplying - H)) 

d 3 q idg f d 3 g 

after an integration by parts over q 1 . The remaining integrals can be solved still 
by integrating by parts over q 1 . The integral proportional to $ ' in Eq. ( 17.71 ) gives 

" q E = - Vi n e , (7.15) 



(2tt) 3 

where we have used the fact that dE/dq l = q l /E. For the integral 

d 3 q q l q k dg 
(2tt) 3 E dq l ' 



(7.16) 



the integration by parts brings two pieces, one from the derivation of q % q k and 
one from the derivation of the energy E 

d 3 q q k f d 3 q q 2 q k k f d 3 q q 2 q k 

/ 7K~^9^^ = -4v n e + / ——g—— .(7.17) 



(27r) 3y E J {2tt) 3 J E E e J (2ir) 3i> E 2 E 

The last integral in Eq. ( |7.17t can indeed be neglected. To check this one makes 
use of the explicit expression (15.41 i for the distribution function g to derive 

dg qi m e 

d^ =9 T e -^^ (7 ' 18) 



and 

d 3 q 



l q j = 5 lJ n e m e T e + n e m 2 e v l v 3 . (7.19) 



(Jhr)3" OT q ' 
Thus it is easy to compute 

* k f d 3 q , 9 T e T e 9 

— Y / TTT^gq \ = -f,kV 2 — + W. k V k n e — + V,kV k V 2 , (7.20) 

which is negligible taking into account that T e /m e is of the order of the thermal 
velocity squared. 

With these results we are now able to compute the left-hand side of the Boltz- 
mann equation (17.8b integrated over d 3 q/(2n) 3 . The same operation must be 
done for the collision terms on the right hand side. For example for the first 
of the equations in (17.8b this brings to the integrals (c ep )QQ> qq i + {c ei ) pp i qq t. 
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However looking at Eq. ( 15.31 ) one realizes that {c ei ) pp i qq > vanishes because the 
integrand is antisymmetric under the change q <-> q' and p <-> p'. In fact this is 
simply a consequence of the fact that the electron number is conserved for this 
process. The same argument holds for the other term (c ep ) qqi qq i . Therefore the 
right-hand side of Eq. ( 17.81 ) integrated over d 3 q/(2n) 3 vanishes and we can give 
the evolution equation for n e . Collecting the results of Eq. ( 17.141 ) to ( 17.201 ) we 
find 

dn^ + ^ + _ + e ^ v k ne ($ _ m k) = Q (7 21) 

OTj OX 

Similarly, for CDM particles, we find 

— q~ — + e Q~i h - * )«CDM 

+ e*+*w£ DM n C DM (V fc - 2* fc ) =0. (7.22) 

7.2. Momentum continuity equations 

Let us now multiply Eq. (|7.7| i by (q l / E) / (2tt) 3 and integrate over d 3 q. In this 
way we will find the continuity equation for the momentum of baryons. The first 
term just gives (n e v 1 )'. The second integral is of the type 

_d_ f totfrf - — (n ^ + n v*iA (7 23) 

dxi J (2n) 39 E E ~ dxi \ e m e + e ) ' ( ' 

where we have used Eq. (17.19b and E = m e . The third term proportional to 

(H - is 

d 3 q u °9 1 l f d 3 q q 2 q l 

f^Se^ = 4n e + / 7Kzb9ipT^ > (7-24) 



(27r) 3 ^ dq k E e J (27r) 3y E 2 E 

where we have integrated by parts over q l . Notice that the last term in Eq. d7.24| i 
is negligible being the same integral we discussewd above in Eq. ( 17.20b . By the 
same arguments that lead to neglect the term of Eq. ( 17.201 ) it is easy to check 
that all the remaining integrals proportional to the gravitational potentials are 
negligible except for 



~ e / 7T^3 ^ = "ee^**' 1 • (7.25) 

J (2ir) s dq k 

The integrals proportional to the second-order vector and tensor perturbations 
vanish as vector and tensor perturbations are traceless and divergence-free. The 
only one which survives is the term proportional to uS 11 + Jiuj 1 in Eq. (17.71 1. 
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Therefore for the integral over d 3 qq l /E of the left-hand side of the Boltz- 
mann equation ( 17.71 ) for a massive particle with mass m e (m p ) and distribution 
function ( 15.41 ) we find 



m e J oxi or] 



(7.26) 



Now, in order to derive the momentum conservation equation for baryons, we 
take the first moment of both Eq. ( 17.81 ) and ( 17.91 ) multiplying them by q and Q 
respectively and integrating over the momenta. Since previously we integrated 
the left-hand side of these equations over d 3 qq l /E, we just need to multiply the 
previous integrals by m e for the electrons and for m p for the protons. Therefore 
if we sum the first moment of Eqs. ( 17.81 ) and ( 17.91 ) the dominant contribution on 
the left-hand side will be that of the protons 

/ T^f® 1 ~dri = lyCepiyql + Q l ^ QQ ' qq ' + ( Ce ^pp'iQ' ■ ( 7 - 27 ) 

Notice that the integral of the Coulomb collision term c ev {q l + Q" 1 ) over all mo- 
menta vanishes simply because of momentum conservation (due to the Dirac 
function 5 4 (q + Q — q' — Q')). As far as the Compton scattering is concerned we 
have that, following Ref. [28], 

{c e ^q )pp r qq' — (c e ^p )pp'qq' , (7.28) 

still because of the total momentum conservation. Therefore what we can com- 
pute now is the integral over all momenta of c ei p l . Notice however that this is 
equivalent just to multiply the Compton collision term C (p) of Eq. (15.31 ) by p 1 
and integrate over d 3 p/ (2n 3 ) 

(ce 7 p l }p P 'qq' = ae* / £Pp*C(p) . (7.29) 



(27T) 

where C(p) has been already computed in Eqs. ( 15.39b and ( 15.40b . 

We will do the integral ( 17.291 ) in the following. First let us introduce the defi- 
nition of the velocity of photons in terms of the distribution function 

(p 7 +p>; = J -^fp\ (7.30) 
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where p 1 = p 7 /3 is the photon pressure and p 1 the energy density. At first order 
we get 

^ = /^A«nS (7.31) 

where A is the photon distribution anisotropies defined in Eq. d6.9l l. At second 
order we instead find 



A ( 2 )« 

4 v\ 
3~ 



i/jA»»*-^. (7.32) 

Therefore the terms in Eqs. ( 15.391 ) and ( 15.40b proportional to / W (p) and /( 2 ) (p) 
will give rise to terms containing the velocity of the photons. On the other hand 
the terms proportional to ffj^ (p) and fffi (p), once integrated, vanish because of 
the integral over the momentum direction n l , J (Kin 1 — 0. Also the integrals 
involving P 2 (v ■ n) = [3(v • n) 2 - l]/2 in the first line of Eq. $5391 and (|5T40T > 
vanish since 

J dOP 2 (v ■ n) n l = v k v j J dflnkUjU 1 = , (7.33) 

where we are using the relations ( 17.13b . Similarly all the terms proportional to 
v, (v • n) 2 and v 2 do not give any contribution to Eq. jl.29i and, in the second- 
order collision term, one can check that j dflY2(ii)n l = 0. Then there are terms 

proportional to (v • n)/(°) (p), (v • n)pdf^ /dp and (v • n)pdf^ / dp for which 
we can use the rules ( I6.10l l when integrating over p while the integration over the 
momentum direction is 



dn, , , r dn 



. v • n)n l = v k / — n K n l = -v l . (7.34) 
47T J Ait 3 

Finally from the second line of Eq. (15.401 1 we get three integrals. One is 

/ p z ( v ' n )/ (1) (p) = h J § A(1) ( v ' n K . ( 7 - 35 > 

where p 1 is the background energy density of the photons. The second comes 
from 



i d& i i k i f d& , i 
Air 3 / 4tt 
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where we have used the rules (I6.10l l. Eq. (17.13b and J (dil/An) n 1, n J n k n l 
{8 ll 8 kl + S lk 5 l: > + d ll S^ k )/l5. In fact the third integral 



exactly cancels the previous one. Summing the various integrals we find 



(7.37) 



dp 

4v (2) 



C(p)p = n e a T p 7 



;(v 



(i) 



a 



iX_ + ^( 1)(v ( 1 )_ v W ) + | *A(')(v.„)n+A«v 



.(7.38) 



Eq. ( 17.38b can be further simplified. Recalling that 5y — A Q we use Eq. ( 17.32b 
and notice that 



where the photon quadrupole IE* is defined as 



(7.39) 



(7.40) 



Thus, our final expression for the integrated collision term ( 17.29b reads 



g(« C1)< -« 7 



4 / «( 2 )* 



„(2)i S 



(7.41) 



We are now able to give the momentum continuity equation for baryons by 
combining m p dg p / drj from Eq. ( 17.26b with the collision term ( 17.291 ) 



d{pbv l ) 
drj 



4(H - ^')p b v r + $' l e 



p b + e + jO b 



T 6 



8 



-n e GTa p-y 



(pbv j v l 



duj 1 
drj 



p b + Huj l pb 



(2)i 
7 



(7.42) 
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where pj, is the baryon energy density and, as we previously explained, we took 
into account that to a good approximation the electrons do not contribute to the 
mass of baryons. In the following we will expand explicitly at first and second- 
order Eq. dT42b . 

7.2. 7. First-order momentum continuity equation for baryons 
At first order we find 



7.2.2. Second-order momentum continuity equation for baryons 
At second order there are various simplifications. In particular notice that the 
term on the right-hand side of Eq. ( 17.42b which is proportional to 6b vanishes 
when matched to expansion of the left-hand side by virtue of the first-order equa- 
tion (17.43b . Thus, at the end we find a very simple equation 




(7.43) 





with t' = — n e ora. 



7.2.3. First-order momentum continuity equation for CDM 
Since CDM particles are collisionless, at first order we find 




(7.45) 



7.2.4. Second-order momentum continuity equation for CDM 
At second order we find 





(7.46) 
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8. Linear solution of the Boltzmann equations 

In this section we will solve the Boltzmann equations at first order in perturbation 
theory. The interested reader will find the extension of these formulae to second 
order in Ref . [2] . The first two moments of the photon Boltzmann equation are 
obtained by integrating Eq. d6.61 > over dfl n /4iv and dfl n n l /4tt respectively and 
they lead to the density and velocity continuity equations 

A oo' + ^v^ l -^ {lY =0, (8.1) 
flyUM* + ±A$'' + = -r 1 (y™ - «W) , (8.2) 



f ^ 4 



where IP 3 is the photon quadrupole moment, defined in Eq. ( 17.40b . 

Let us recall here that = A$ = / df2A (1) /47r and that the photon 
velocity is given by Eq. (17.311 1. 

The two equations above are complemented by the momentum continuity 
equation for baryons, which can be conveniently written as 



T 



(8.3) 



where we have introduced the baryon-photon ratio R = 3pb/(4p 7 ). 

Eq. (18.31 ) is in a form ready for a consistent expansion in the small quantity r~ 1 



which can be performed in the tight-coupling limit. By first taking ijW* — v^ 1 
at zero order and then using this relation in the left-hand side of Eq. (18.31 ) one 
obtains 



7 T 



(8.4) 



Such an expression for the difference of velocities can be used in Eq. ( 18.21 ) to give 
the evolution equation for the photon velocity in the limit of tight coupling 

+ u iTr v( ' )i + \t+~r + = ■ (8 ' 5) 



Notice that in Eq. ( 18.51 ) we are neglecting the quadrupole of the photon distri- 
bution n^ 1 ^ (and all the higher moments) since it is well known that at linear 
order such moment(s) are suppressed in the tight-coupling limit by (successive 
powers of) 1 /r with respect to the first two moments, the photon energy den- 
sity and velocity. Eqs. ( 18. U and ( 18.5b are the master equations which govern 
the photon-baryon fluid acoustic oscillations before the epoch of recombination 
when photons and baryons are tightly coupled by Compton scattering. 
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In fact one can combine these two equations to get a single second-order dif- 
renti; 

yields 



ferential equation for the photon energy density perturbations A^. Deriving 
Eq. (18. lb with respect to conformal time and using Eq. (18.5b to replace diVj 



-c 2 V 2 (Ag> - 4*«) = ^V 2 + ^) , (8.6) 



where c s = 1/ -^/3(1 + i?) is the speed of sound of the photon-baryon fluid. 
Indeed, in order to solve Eq. (18.6b one needs to know the evolution of the grav- 
itational potentials. We will come back later to the discussion of the solution of 
Eq. (US). 

A useful relation we will use in the following is obtained by considering the 
continuity equation for the baryon density perturbation. By perturbing at first 
order Eq. ( 17.21) we obtain 

<5p + v*i - = . (8.7) 

Subtracting Eq. ( 18.7b form Eq. ( 18.1b brings 

A oo'-^ 1) '+ ^-v^) ti = 0, (8.8) 
which implies that at lowest order in the tight-coupling approximation 

A« = ~/ b 1] , (8-9) 
for adiabatic perturbations. 

8.1. Linear solutions in the limit of tight coupling 

In this section we briefly recall how to obtain at linear order the solutions of 
the Boltzmann equations ( 18.6b . These correspond to the acoustic oscillations of 
the photon-baryon fluid for modes which are within the horizon at the time of 
recombination. It is well known that, in the variable (Ago — 4\I/( 1 )), the solution 
can be written as [32, 38] 

[1 + i?(?y)] 1 / 4 (A[ ) 1 ) - 4*< 1 )) = Acos[kr s (r))} + B sin[fcr s (r?)] 



-4- 



XfflnMr.faJ-r.fa 7 ))], (8.10) 
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where the sound horizon is given by r s (77) = J? drj'c s (i]'), with R = 3/9(,/(4p 7 ). 
The constants A and B in Eq. (18.10b are fixed by the choice of initial conditions. 

In order to give an analytical solution that catches most of the physics under- 
lying Eq. (18.10b and which remains at the same time very simple to treat, we will 
make some simplifications following Ref. [28,39]. First, for simplicity, we are 
going to neglect the ratio R wherever it appears, except in the arguments of the 
varying cosines and sines, where we will treat R — R^asa constant evaluated at 
the time of recombination. In this way we keep track of a damping of the photon 
velocity amplitude with respect to the case R = which prevents the acoustic 
peaks in the power-spectrum to disappear. Treating R as a constant is justified 
by the fact that for modes within the horizon the time scale of the ocillations is 
much shorter than the time scale on which R varies. If R is a constant the sound 
speed is just a constant c s = 1/^3(1 + i?*), and the sound horizon is simply 
r s {ri) = c s r\. 

Second, we are going to solve for the evolutions of the perturbations in two 
well distinguished limiting regimes. One regime is for those perturbations which 
enter the Hubble radius when matter is the dominant component, that is at times 
much bigger than the equality epoch, with k <C k eq ~ Tleq> where k eq is the 
wavenumber of the Hubble radius at the equality epoch. The other regime is 
for those perturbations with much smaller wavelenghts which enter the Hubble 
radius when the universe is still radiation dominated, that is perturbations with 
wavenumbers k ^> k eq ~ rj^ 1 . In fact we are interested in perturbation modes 
which are within the horizon by the time of recombination 77* . Therefore we will 
further suppose that 77* S> rj eq in order to study such modes in the first regime. 
Even though 77* ^> r) eq is not the real case, it allows to obtain some analytical 
expressions. 

Before solving for these two regimes let us fix our initial conditions, which are 
taken on large scales deep in the radiation dominated era (for r\ — ► 0). During this 
epoch, for adiabatic perturbations, the gravitational potentials remain constant on 
large scales (we are neglecting anisotropic stresses so that "J)' 1 ) ~ 1 J?( 1 )) and from 
the (0 — 0)-component of Einstein equations 



On the other hand, from the energy continuity equation (18.1b on large scales 



from Eq. ( 18.1 lb the constant on the right-hand side of Eq. ( 18.121 ) is fixed to be 
-6$W(0); thus we find B = and A = -6$W(0). 



$«(0)=-±AW(0). 



(8.11) 



Ajjo 3 - 4* (1) = const. ; 



(8.12) 
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With our simplifications Eq. ( 18.10b then reads 

A^ - 4* (1) = - 6$ (1) (0)cos(w 7?) 

- drf^ l \rf)an[uo(v-rf)], (8.13) 

where uq — kc s . 

8.2. Perturbation modes with k <C k eq 

This regime corresponds to perturbation modes which enter the Hubble radius 
when the universe is matter dominated at times rj ^> rj eq . During matter dom- 
ination the gravitational potential remains constant (both on super-horizon and 
sub-horizon scales), as one can see for example from Eq. ( IB. Il l, and its value 
is fixed to $^(k, rj) = ^$^(0), where $W(0) corresponds to the gravita- 
tional potential on large scales during the radiation dominated epoch. Since we 
are interested in the photon anisotropics around the time of recombination, when 
matter is dominating, we can perform the integral appearing in Eq. ( 18.101 ) by tak- 
ing the gravitational potential equal to its value during matter domination so that 
it is easily computed 

r v 18$ (1) (0) 

2 / dr/SW^sinNfr-r/)] = ° W (i _ cos (cu aV )) . (8.14) 

Jo 10 ujq 

Thus Eq. ( 18. 1 3b gives 

A$ -4*^ = 5$(i)(o)cos(a;o»7) - ^$ (1) (0) . (8.15) 

The baryon-photon fluid velocity can then be obtained as diV^ 1 = — 3(A^^ 
4^r(i))'/4 f rom Eq. (|QT) , In Fourier space 



on 



ih = ^-$ (1) (0) sin^H , (8.16) 



10 



where, going to Fourier space, diV$p % — » iki uiy (k) and 

v i )l = -ij^ {1) (0)sm(uQV)c a , (8.17) 



since the linear velocity is irrotational. 
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8.3. Perturbation modes with k 3> k eq 

This regime corresponds to perturbation modes which enter the Hubble radius 
when the universe is still radiation dominated at times 77 <C rj eq . In this case 
an approximate analytical solution for the evolution of the perturbations can be 
obtained by considering the gravitational potential for a pure radiation dominated 
epoch, given by Eq. ( IB. 81 ). For the integral in Eq. ( 18.131 ) we thus find 

[\^(r]')sm[uj (ri - r, 1 )] = - JL cos ( Wo7? ) , (8.18) 

where we have kept only the dominant contribution oscillating in time, while 
neglecting terms which decay in time. The solution d8.13t becomes 

A$ - 4* (1) = 6$ (1) (0)cos(w ry) , (8.19) 
and the velocity is given by 

v -' )i = ~4^ $(1)(o) sin ^ cs > (8 - 2o) 

Notice that the solutions (18.1 9b — (18 . 20b are actually correct only when radiation 
dominates. Indeed, between the epoch of equality and recombination, matter 
starts to dominate. Full account of such a period is given e.g. in Section 7.3 of 
Ref. [28], while its consequences for the CMB anisotropy evolution can be found 
e.g. inRef. [40]. 



9. Conclusions 

In these lecture notes we derived the equations which allow to evaluate CMB 
anisotropies, by computing the Boltzmann equations describing the evolution of 
the baryon-photon fluid up to second order. This allows to follow the time evo- 
lution of CMB anisotropies (up to second order) on all scales, from the early 
epoch, when the cosmological perturbations were generated, to the present time, 
through the recombination era. The dynamics at second order is particularly im- 
portant when dealing with the issue of non-Gaussianity in CMB anisotropies. 
Indeed, many mechanisms for the generation of the primordial inhomogeneities 
predict a level of non-Gaussianity in the curvature perturbation which might be 
detectable by present and future experiments. 
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Symbol Definition Equation 



<E>, ^ Gravitational potentials in Poisson gauge 13. U 

u>i 2nd-order vector perturbation in Poisson gauge I3.lt 

Xij 2nd-order tensor perturbation in Poisson gauge I3.lt 

r) Conformal time 13. U 

/ Photon distribution function 14. 81 

q Distribution function for massive particles H3} & ED 

fW i-th order perturbation of the photon distribution function 14.231 

Moments of the photon distribution function )5. 171 

C(p) Collision term (33) & (5J\ 

p Magnitude of photon momentum (p = pn') 14.3) 

n 1 Propagation direction 14.6) 

A' 1 ) (x l , n' , rf) First-order fractional energy photon fluctuations <6.5t 

A"' (x l , n' , rj) Second-order fractional energy photon fluctuations 46.9t 

n e Electron number density 17.1 li 

S e (5f,) Electron (baryon) density perturbation 15.131 

k Wavenumber 16.16) 

v m Baryon velocity perturbation d6.19) & 16.20) 

^CDM Cold dark matter velocity 17.46) 

Second-order photon velocity {1. 321 

Si m Temperature source term 16.17) 

r Optical depth 16.14) 

p 7 (pi,) Background photon (baryon) energy density 17.43) 



Appendix A. Einstein's equations 

In this Appendix we provide the necessary expressions to deal with the gravita- 
tional part of the problem we are interested in, namely the second-order CMB 
anisotropics generated at recombination as well as the acoustic oscillations of 
the baryon-photon fluid. The first part of the Appendix contains the expressions 
for the metric, connection coefficients and Einstein's tensor perturbed up to sec- 
ond order around a flat Friedmann-Robertson- Walker background, the energy- 
momentum tensors for massless (photons) and massive (baryons and cold dark 
matter) particles, and the relevant Einstein's equations. The second part deals 
with the evolution equations and the solutions for the second-order gravitational 
potentials in the Poisson gauge. 

Appendix A. 1. The metric tensor 

As discussed in Section 3, we write our second-order metric in the Poisson gauge, 

ds 2 = a 2 (r]) [-e^df] 2 + 2uj t dx l dr] + (e~ 2 %- + Xl] )dx l dx j } , (A.l) 
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where 0(77) is the scale factor as a function of conformal time 77, and Ui and Xij 
are vector and tensor peturbation modes respectively. Each metric perturbation 
is expanded into a linear (first-order) and a second-order part, as discussed in 
Section 3. 

Appendix A.2. The connection coefficients 
For the connection coefficients we find 



1 00 

1 Oi 



Oj 



jk 



= H + &, 



1 00 — 



+ 



2*+2$ 



dxi 



1 / <9w, 



2 
a* 



1 / d x \ | d x \ d Xj k 

2 \ 3x fe dxi dxi 



(A.2) 



Appendix A. 3. Einstein tensor 

The components of Einstein's tensor read 

-2* 



G° n = -- 



+ 



G\ = 



2$+2* 



[m 2 - 6H^' + 3(tf ') 
(d^d^ - 2V 2 *)] , 



1 



-2$ 



+ 2f" 



2^ 


)- 


a , 


' a 2 ' 








a" 

_2— -2*'$' 




a 


e 2 *( 


<9 fc $<9 fe $ + V 2 $ 



(A.3) 



(A.4) 



,/\2 



„23< 
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H . „ .-s 1 



(d^ ] +d J ^)- — (dW J + d J J) 



+ ^(^ + lxf-^ 2 xy). (A.5) 

Taking the traceless part of Eq. ( IA.5b . we find VP — $ = Q, where Q is defined 

by V 2 Q = -P + 3N, with P = P l u 

P* = a*$a,-* + - (5**3,-* - a**©,-*) + 47rG N a 2 e- 2 *:r ; - (A.6) 

and V 2 iV = djd j PL 

The trace of Eq. jA.St gives 

e -2$ / ^2 _ 2 — - 2$'*' - 3(#') 2 + (3*' - Q') + 2*" 



e 2 * 



(2<9 fc $<9 fc $ + a fc *<9 fe * - 2<9 fe $d fc * + 2(P - 3JV)) 

3 

= ^« a T*. (A.7) 
From Eq. (1A.4K we may deduce an equation for uo l 

- ^VV + (w 2 - 2^-\ uj 1 (A.8) 



Si - 1^1 

°3 y 2 



(2 (W + (W - *') - 87rG N a 2 e- 2 *T J ) 



Appendix A.4. Energy-momentum tensor 

Appendix A.4.1. Energy-momentum tensor for photons 
The energy-momentum tensor for photons is defined as 



2 r d 3 P p^p u 

= g] (2^)3 ~P°" 



where g is the determinant of the metric ( I A. It and / is the distribution function. 
We thus obtain 

T% = -pJ l + A^ + ^f] , (A.10) 
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T* o = -\^ Pl + + A#„W*) + \- Pl e*-*J (A.ll) 

t; 3 = a, fn; i + irf, (i + a« + H , (A.i2) 

where p 7 is the background energy density of photons and 
is the quadrupole moment of the photons. 

Appendix A.4.2. Energy-momentum tensor for massive particles 

The energy-momentum tensor for massive particles of mass m, number density 

n and degrees of freedom gd 

l ™» " ^Tg-J ( 27r )3 Q0 3m, IA.14J 

where g m is the distribution function. We obtain 

T^ = -^m = -Pm(l + ^ ) + ^ ) ) , ( A -!5) 

= — e 



= -e*+*p m + + (A.16) 

= Pro p~ + i ) = P m I 5 J ■ ~ + U ™ % j 



i = Pm I <5j + <«m * I = Pm I ^ + <^^j , (A. 17) 

where p m is the background energy density of massive particles and we have 
included the equilibrium temperature T m . 



Appendix B. First-order solutions of Einstein's equations in various eras 

Appendix B.l. Matter-dominated era 

During the phase in which the CDM is dominating the energy density of the Uni- 
verse, a ~ ?/ 2 and we may use Eq. ( IA.7b to obtain an equation for the gravitational 
potential at first order in perturbation theory (for which = <3/( 1 )) 

+37i$ (1) ' = 0, (B.l) 
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which has two solutions = constant and = H/a 2 . At the same order 
of perturbation theory, the CDM velocity can be read off from Eq. iAAi 

v Wi = — . (B.2) 
The matter density contrast 5^ satisfies the first-order continuity equation 

Going to Fourier space, this implies that 

where 5^ (0) is the initial condition in the matter-dominated period. 



Appendix B. 2. Radiation-dominated era 

We consider a universe dominated by photons and massless neutrinos. The 
energy-momentum tensor for massless neutrinos has the same form as that for 
photons. During the phase in which radiation is dominating the energy density 
of the Universe, a ~ 77 and we may combine Eqs. ( IA.3I ) and ( IA.7I ) to obtain an 
equation for the gravitational potential at first order in perturbation theory 



4W *<i)'_Iv 2 *< 1 > = WQ (1) ' + iv 2 Q«, 
3 3 



v2g( i) = 9 H 2^ n (i)^ (BS) 
where the total anisotropic stress tensor is 



n^ = ^n; J + ^n^. (B.6) 

PT PT 

We may safely neglect the quadrupole and solve Eq. (IB. 51 setting u± = 



Then Eq. ( |B.5t , in Fourier space, becomes 
2 

3 fj* 



«" + -«'+ f^- 4 l« = 0. (B.7) 
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This equation has as independent solutions u + = ji(kr)/\/S), the spherical 
Bessel function of order 1, and w_ = n\{krj/\/i), the spherical Neumann func- 
tion of order 1 . The latter blows up as r\ gets small and we discard it on the basis 
of initial conditions. The final solution is therefore 

T m „ t(i1 . ,sm 
k ~ [ 1 (k V /V3) 3 

where $W (0) represents the initial condition deep in the radiation era. 

At the same order in perturbation theory, the radiation velocity can be read off 
from Eq. ( lA.4b 

< } = " — L ■ (B.9) 



References 

[1] N. Bartolo, S. Matarrese and A. Riotto, JCAP 0606, 024 (2006) |arXiv:astro-ph/0604416| . 
[2] N. Bartolo, S. Matarrese and A. Riotto, JCAP 0701, 019 (2007) [arXiv:astro-ph/0610110]. 
[3] For a review, see D. H. Lyth and A. Riotto, Phys. Rept. 314, 1 (1999). 
[4] D. N. Spergel et at , |arXiv:astr7>p h/0603449 
[5] See http://plcmck.esa.int/ 

[6] N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, Phys. Rept. 402, 103 (2004) 

[arXiv:astro-ph/0406398~| . 
[7] V. Acquaviva, N. Bartolo, S. Matarrese and A. Riotto, Nucl. Phys. B 667, 1 19 (2003); J. Mal- 

dacena, JHEP 0305, 013 (2003). 
[8] N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. D 65, 103505 (2002); F. Bernardeau and 

J. P. Uzan, Phys. Rev. D 66, 103506 (2002); F. Vernizzi and D. Wands, |arXTv:astro- ph/0603799 
[9] D. H. Lyth, C. Ungarelli and D. Wands, Phys. Rev. D 67 (2003) 023503. 
[10] T. Hamazaki and H. Kodama, Prog. Theor. Phys. 96 (1996) 1 123. 

[11] G. Dvali, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 69, 023505 (2004); L. Kofman, 

|arXiv:astro-ph/03036T4l 
[12] G. Dvali, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 69 (2004) 083505. 
[13] C. W. Bauer, M. L. Graesser and M. P. Salem, Phys. Rev. D 72 (2005) 023512. 
[14] D. H. Lyth, JCAP 0511 (2005) 006; M. P. Salem, Phys. Rev. D 72 (2005) 123516. 
[15] M. Bastero-Gil, V. Di Clemente and S. F. King, Phys. Rev. D 70, 023501 (2004); E. W. Kolb, 

A. Riotto and A. Vallinotto, Phys. Rev. D 71, 0435 13 (2005); C. T. Byrnes and D. Wands, Phys. 

Rev. D 73, 063509 (2006) |arXiv:astro-ph/0512195| . 
[16] T. Matsuda, Phys. Rev. D 72 (2005) 123508. 

[17] W. Hu, Phys. Rev. D 64, 083005 (2001); T. Okamoto and W. Hu, Phys. Rev. D 66, 063008 
(2002). 

[18] G. De Troia et at, Mon. Not. Roy. Astron. Soc. 343, 284 (2003). 

[19] J. Lesgourgues, M. Liguori, S. Matarrese and A. Riotto, Phys. Rev. D 71, 103514 (2005). 



Cosmic Microwave Background Anisotropics up to Second Order 



53 



[20] N. Bartolo, S. Matarrese and A. Riotto, JCAP 0605, 010 (2006) |arXiv:astro-ph/05 12481 |. 

[21] T. Pyne and S. M. Carroll, Phys. Rev. D S3, 2920 (1996); S. Mollerach and S. Matarrese, Phys. 
Rev. D 56, 4494 (1997); S. Matarrese, S. Mollerach and M. Bruni, Phys. Rev. D 58, 043504 
(1998); N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. Lett. 93, 231301 (2004); K. Tomita, 
Phys. Rev. D 71, 083504 (2005); N. Bartolo, S. Matarrese and A. Riotto, JCAP 0508, 010 
(2005). 

[22] W. Hu and S. Dodelson, "Cosmic Microwave Background Anisotropics, Ann. Rev. Astron. 
Astrophys. 40, 171 (2002). 

[23] U. Seljak and M. Zaldarriaga, "Direct Signature of Evolving Gravitational Potential from Cos- 
mic Microwave Background," Phys. Rev. D 60, 043504 (1999). 

[24] D. M. Goldberg and D. N. Spergel, "Microwave Background Bispectrum. 2. A Probe Of The 
Low Redshift Universe," Phys. Rev. D 59, 103002 (1999). 

[25] D. M. Goldberg and D. N. Spergel, Phys. Rev. D 59, 103002 (1999). 

[26] E. Komatsu and D. N. Spergel, Phys. Rev. D 63, 063002 (2001) 

[27] P. Creminelli and M. Zaldarriaga, Phys. Rev. D 70, 083532 (2004). 

[28] S. Dodelson, Modern cosmology, Amsterdam, Netherlands: Academic Pr. (2003) 440 pp. 

[29] S. Matarrese, S. Mollerach and M. Bruni, Phys. Rev. D 58, 043504 (1998). 

[30] W. Hu, D. Scott and J. Silk, Phys. Rev. D 49, 648 (1994). 

[31] S. Dodelson and J. M. Jubas, Astrophys. J. 439, 503 (1995). 

[32] W. T. Hu, PhD thesis Wandering in the background: A Cosmic microwave background explorer, 

|arXiv;astro- ph/9508 1 26 
[33] R. Maartens, T. Gebbie and GR. Ellis, Phys. Rev. D 59, 083506 (1999). 
[34] W. Hu, U. Seljak, M. J. White and M. Zaldarriaga, Phys. Rev. D 57, 3290 (1998) 
[35] U. Seljak and M. Zaldarriaga, Astrophys. J. 469, 437 (1996) 
[36] W. Hu and M. J. White, Phys. Rev. D 56, 596 (1997) 
[37] W. Hu, Astrophys. J. 529, 12 (2000) 
[38] W. Hu and N. Sugiyama, Astrophys. J. 444, 489 (1995) 
[39] M. Zaldarriaga and D. D. Harari, Phys. Rev. D 52, 3276 (1995) 
[40] V. F. Mukhanov, Int. J. Theor. Phys. 43, 623 (2004) 



